Evaluating individual variability in species‐habitat relationships using an integrated step-selection analysis
We provide worked-out examples, code templates, and functions for quantifying variability in movement and habitat-selection parameters when implementing Integrated Step-Selection Analysis with random effects.
Integrated step-selection analyses (ISSAs) (Avgar et al. 2016; Fieberg et al. 2021) are frequently used to study animal movement and habitat selection. Methods for incorporating random effects in ISSAs have been developed, making it possible to quantify variability among animals in their space-use patterns. Although it is possible to model variability in both habitat selection and movement parameters (Muff, Signer, and Fieberg 2020), applications to date have focused on the former despite the widely acknowledged and important role that movement plays in determining ecological processes from the individual to ecosystem level. One potential explanation for this omission is the absence of readily available software for estimating movement parameters in ISSAs with random effects. We demonstrate methods for estimating individual-specific movement parameters, including situations where step-lengths and turn angles depend on continuous or categorical habitat variables.
\[\begin{equation} u(s, t+ \Delta t)|u(s',t), \psi = \frac{w(x(s,t+\Delta t); \beta)\phi(s,s', x(s', t);\gamma)}{\int_{\tilde{s}\in\ G} w(x(\tilde{s}, t+\Delta t); \beta)\phi(\tilde{s}, s', x(s', t);\gamma)d\tilde{s}}, \tag{1} \end{equation}\]
Movement-free selection kernel, (\(w(.)\)), was modeled as a log-linear function of reef class and distance-to-edge measured at the end of the movement step.
We fit two different habitat-dependent movement kernels,
The number of relocations (range 265-52,189) as well as the number of tracking days (range 3-223) were substantially different across individuals (n=35). Coefficients from the integrated step-selection models were used to understand how individuals move through different types of habitat using the selection-free movement kernels.There was substantial variability in the estimated individual-specific habitat selection parameters, with both positive and negative coefficients for each explanatory variable (Figure 3). We plotted the step-length and turn-angle distribution for all the individuals (blue) and a “typical individual”, i.e., one with all random effects equal to 0 (Figure 1, Figure 2). Snappers move shorter distances, on average, when in high-relief habitat and close to reef edges. Also, their movements were least directed (with more dispersed turn angles) when they were in sand or low-relief habitats or away from the reef edges .
The typical individual strongly selects for high-relief hardbottom habitat relative to sandy habitat and for areas close to the edge of sand/low-relief and medium/high-relief hardbottom classes (Figure 3). Individuals also prefer areas close to the reef.
Figure 1: Estimated step-length distributions for different individuals and for a typical individual in different habitats (A-D). Blue lines show the step-length distribution for individual fish and red line shows the distribution for a typical indiviudal (i.e., an individual with all random effects set to 0).
Figure 2: Estimated turn-angle distributions for different individuals and for a typical individual in different habitats (A-D). Blue lines show the step-length distribution for individual fish and red line shows the distribution for a typical indiviudal (i.e., an individual with all random effects set to 0).
Figure 3: Boxplot summarizing the individual habitat-selection coefficients (black dots) associated with reef class indicator variables and the distance-to-edge covariate for Model 1 (A) and Model 2 (B). Boxes bound the 25th and 75th percentiles, solid line within the box indicates the median, and the whiskers extend to 1.5 times the interquartile range of the observations. The mean coefficient, along with a 95% confidence interval for the mean, is depicted in brown.